Heat flow in InAs/InP heterostructure nanowires 
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The transfer of heat between electrons and phonons plays a key role for thermal management in 
future nanowire-based devices, but only a few experimental measurements of electron-phonon (e-ph) 
coupling in nanowires are available. Here, we combine experimental temperature measurements on 
an InAs/InP heterostructure nanowire system with finite element modeling (FEM) to extract infor- 
mation on heat flow mediated by e-ph coupling. We find that the electron and phonon temperatures 
in our system are highly coupled even at temperatures as low as 2 K. Additionally, we find evidence 
that the usual power-law temperature dependence of electron-phonon coupling may not correctly 
describe the coupling in nanowires and show that this result is consistent with previous research on 
similar one-dimensional electron systems. We also compare the strength of the observed e-ph cou- 
pling to a theoretical analysis of e-ph interaction in InAs nanowires, which predicts a significantly 
weaker coupling strength than observed experimentally. 

PACS numbers: 44.10.-K, 63.20.kd, 66.30.Pa 



I. INTRODUCTION 

There exists great potential for the use of nanowires 
(NW) for future nanoelectronic applications, such 
as light-emitting diodes (LEDs)ii2, photovoltaic (PV) 
cells 3 -^, wrap-gate transistors 5 -^, and low-dimensional 
thermoelectric^—. In many of these devices heat flow 
plays a significant role in device performance, either be- 
cause heat flow is a parasitic effect and is therefore un- 
desired, as is the case in thermoelectrics, or because 
high heat flow is critically required for thermal manage- 
ment, as in most other applications, such as LEDs, PV 
cells and transistors. In many cases, heat is produced 
in the form of Joule heat in the electronic system, but 
is then distributed between the electrons and phonons 
through electron-phonon coupling. Since phonons gener- 
ally have the higher thermal conductivity in semiconduc- 
tors, electron-phonon coupling has a central role in heat 
flow through semiconductor nanowires. 

Few experimental papers discuss e-ph interaction in 
semiconductor nanowire systems, where electrons and 
phonons are confined to one-dimensional transport at 
low temperatures On the theory side, e-ph inter- 
action strength was previously analyzed using three- 
dimensional (3D) phonon systems, which typically leads 
to the prediction of T 3 or T 5 power laws depending 
on the type of coupling considered] 14 ' 15 However, in ex- 
perimental systems with ID electrons and 3D phonons, 
experimentally observed power laws do not agree with 
these predictions J^r— In vertically grown NWs then, 
where both electrons and phonons can be confined to 
ID, one would expect the same, if not greater, devia- 
tion from the usual power laws. In fact, many previous 
theoretical analyses of phonon confinement in nanoscale 
structures have predicted increased e-ph interaction en- 



ergy exchange and scattering rates compared to bulk 
systems^—, deviations from the T 3 or T 5 power laws 2 ^, 
and e-ph interaction regimes where powers laws are not 
obeyed^ 3 .. Together, these previous results highlight 
the need for additional investigation of e-ph coupling in 
nanowires, both experimentally and from a fundamental 
point of view. 



To probe e-ph coupling in NWs, we have carried out 
temperature measurements of the electron temperature 
profile around a double-barrier heterostructure embed- 
ded into an InAs nanowire (Fig. [T]). When one end of the 
nanowire is heated, we observe an increase in the elec- 
tron temperature on the cold side of the nanowire, which 
is not expected when considering only electronic diffu- 
sive heat flow. By combining these measurements with a 
FEM model of our heterostructure nanowire (HNW) sys- 
tem, we extract information about e-ph coupling, as well 
as other thermal transport properties, within the InAs 
portion of the HNW. Specifically, we find that e-ph in- 
teraction within the HNW plays a key role in explaining 
the observed electron temperatures, and that we must 
consider both the electronic and phononic thermal con- 
ductivities. 



In the following, we first describe our heterostructure 
nanowire system and temperature measurements. We 
follow this with a discussion on our heat flow model for 
the HNW and its implementation using FEM modeling. 
Following a discussion of the FEM results, we lay out 
theoretical calculations of e-ph interaction in InAs NWs. 
Lastly, we compare our FEM and theory results on e-ph 
interaction in NWs with literature values. 
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FIG. 1. (a) A scanning electron micrograph (SEM) image of 
the HNW system showing the hot (HC) and cold (CC) Au 
contacts, and the approximate position of the quantum dot 
(QD), not visible here. A temperature gradient is applied us- 
ing an AC heating current, In, through the HC. (b) Close-up 
SEM image of the two InP barriers within the InAs nanowire. 
(c) Cross section of the nanowire and contacts along the 
nanowire axis. AT( H c,cc) are the heating-induced electron 
temperature rises above the background temperature, To, in 
the hot and cold contacts respectively. AT ej ( Si( j) are the source 
and drain electron temperature rises measured in the vicin- 
ity of the quantum dot. The source and drain portions of the 
NW are defined by the direction of heat flow, the source being 
closest to the HC. 



II. EXPERIMENTAL DEVICE AND 
TEMPERATURE MEASUREMENTS 

The HNW contains a quantum dot (QD) defined by 
two 4 nm InP barriers, spaced approximately 14 nm 
apart (Fig. [2(b)). The HNW has an average diameter 
of 55 nm, a length of 1.26 /im between its two Ni/Au 
(25/75 nm) contacts, and an estimated carrier density of 
n C arrier = 2.6 ± 0.3 x 10 17 cm -3 at 4.2 K. The leads and 
HNW lie on top of a Si substrate with a 100 nm thick 
SiC>2 capping layer. The measurement set-up and ther- 
mometry techniques are described in detail in Refs. [24] 
and [25]. In brief, two 180° out-of-phase AC heating volt- 
ages applied to the hot contact (HC) are used to heat 
the electrons on the source side of the HNW (Fig. [T|). 
Note that the heating current (Jr) does not actually flow 
through the HNW itself. This is achieved by tuning the 
amplitudes of the two voltages to create an AC voltage 
node at the location of the HNW, thus keeping the bias 
voltage due to the heating current at negligible levels. 
The temperature of electrons entering the QD is deter- 
mined by comparing temperature-bias-driven current to 
voltage-bias-driven current i 24 i 25 An external bias voltage 
applied across the HNW controls the direction of cur- 
rent flow, which determines whether AT c .s or AT e ^ (see 
Fig. [1} is measured. If electrons flow from the HC to the 
CC, AT 6jS is measured, while electron flow from the CC 
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FIG. 2. Measured source (red filled circles), AT CjS , and drain 
(blue open circles), AT Ci d, electron temperature rises for To 
= 1.228, 2.2, 2.94, and 4.25 K. The relatively large tempera- 
ture rise observed for the drain electrons, combined with the 
low electron conductance of the QD, indicates the presence of 
phonon-mediated heat flow into the drain electron reservoir. 



to the HC measures AT eA H& 

In Fig.[5]we show measurements of AT e ( s d ) (In) at four 
different background temperatures, T$: 1.228, 2.2, 2.94, 
and 4.25 K. Our HNW has an electrical resistance of 3.2 
Mfi, roughly three orders of magnitude higher than in 
similar InAs NWs without embedded hetcrostructures. 
We can therefore assume that the electric resistance is 
dominated by the QD. Within an electronic diffusive heat 
model, the QD's high electric resistance would lead one 
to expect virtually all of the temperature differential ap- 
plied to the HNW to fall across the QD. If this model were 
complete, the drain electron temperature should remain 
at the background temperature, AT c c i w 0. In contrast, 
for each To we find an unexpectedly warm drain elec- 
tron temperature, suggesting the presence of at least one 
significant, additional heat flow mechanism warming the 
drain electron reservoir. Through FEM modeling, we will 
show that e-ph coupling within the HNW can explain our 
key observation of an increased AT Ci d- 

III. HNW HEAT FLOWS AND FINITE 
ELEMENT MODELING 

To study heat flows through the HNW system, we 
separately model the HNW and the surrounding leads 
and substrate using finite element modeling in COMSOL 
Multiphysics, see Fig. [3j We model the electron, T , and 
phonon, T p h, temperatures using two coupled diffusive 
heat equations, 

-V • (k c (T c ) VT C ) = Qj - Q c _ ph , (1) 
-V • (ftphCTph) VT ph ) = Qe-ph ■ (2) 

Here K( e ,ph) are the electron and phonon thermal conduc- 
tivities, respectively, and Qj represents Joule heat gen- 
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FIG. 3. (a,b) Top and (c) side view schematics of the simu- 
lated area for the Au leads and substrate, called the substrate 
simulations in the text. The location of where the HNW 
would be is shown by dashed lines in panels (b) and (c). The 
outer boundary in (a) and lower boundaries in (c) (blue lines) 
are set to the background temperature, To- The top boundary 
in (c) (purple line) is assumed to be a perfect insulator. The 
red line in panel (b) shows the path of the heating current, Ih- 
Panel (d) shows a schematic of the simulated HNW. In panels 



(c) and (d) , AT^F and ATq^F represent the electron and 
phonon temperature rises in the HC and CC, respectively, as 
well as the boundary conditions for the simulated HNW. 



crated in the electron system. The Q c - P h term describes 
the heat exchanged between electrons and phonons (e-ph 
interaction) in the Au leads and HNW, and we assume 
it follows a general power law of the form, 



Q 



e-ph 



r(TT 



ph) 



(3) 



where the parameters T and n describe the strength and 
type of e-ph coupling respectively^ 

We simulate the electric potential, 0, through the 
metallic leads using Laplace's equation, 



V^ = 0, 



(4) 



with the primary purpose of determining the amount of 
Joule heating in the hot contact. However, to determine 
the amount of Joule heating in the HNW, we use the 
experimentally measured electric current and electrical 
conductance of the HNW. 

We account for the temperature dependencies of n e and 
ftph by using the Wiedemann-Franz law for electrons, 



KeC^e) = L aT c 



(5) 



and a cubic power law for phonons, valid well below the 
Debye temperature, 



K ph(?ph) 



(6) 



Here, L Q = (n 2 /3)(fc B /e) 2 = 2.44 x 1CT 8 Wfl/K 2 is 
the Lorenz number, a is the electrical conductivity, and 



Parameter 


Au 


C ph (W/K 4 m) 


3 x 10 -3W 


C c (1/ftm) 


4 x 10 7(b) 


T (W/m 3 K 5 ) 


1.8 x 10 9(c) 


n 


5 «3 



a Bulk values^ for Au assuming a phonon mean free path of 

100 nm, the film thickness. 
b Measured separately on a Au thin film structure similar to 

those used here. 
c Average of the values from Refs. [23] andl2gl. 

TABLE I. Simulation parameters used for the Au leads. 



Parameter 


Lower 
Bound 


Upper 
Bound 


Comments 


Cph 
(W/K 4 m) 


5 x 1(T 5 


1.5 x 10~ 3 


See footnote (jsj 


C c 
(1/fim) 


5 x 10 2 


5 x 10 5 


Measurements 
on pure InAs NWs^ 


r 

(W/m 3 K™) 


10 6 


10 13 


Refs. 11, 12, 13, 
27, 29-31 03 


n 


0.1 


7 



a Values initially based on Drude model, but were adjusted to 
conform to a region where the FEM simulations converge. 

b Ranges extended beyond values found in references to 
encompass possible deviations. 

c References include examples of bulk, thin film, and NW values. 

TABLE II. HNW parameter ranges used in the Nelder-Mead 
optimization method. 



Cph is a material constant. We assume that the De- 
bye temperatures can be described by their bulk values: 
e£ u 163 K, and 6^ lAs ps 255 K. In general, a could be 
temperature dependent, 



a = C e f(T e ) 



(7) 



where f(T) is some temperature dependent, dimension- 
less function; however, based on electrical conductivity 
measurements on similar InAs nanowires and thin film 
Au structures, we assume f(T) = 1 for the range of tem- 
peratures looked at here. 

For the surrounding leads to the HNW, we assume that 
the transport properties are well described by bulk Au, 
and will stop mentioning the Ni component. The values 
we use for the four available parameters, {C , C p h, I\ n}, 
for the Au leads are shown in Table HI Unfortunately, the 
setup of our system prevents us from measuring any of 
these four parameters for the InAs portion of the HNW 
directly. However, based on prior research on similar 
InAs NWs and e-ph coupling, we can define a range of 
reasonable values for our HNW, see Table HU Tabulated 
values from Refs. HI and HH are used for the Si and SiC>2 
phonon thermal conductivities. We assume no electric 
current flows through the Si and SiC>2. 

During thermometry measurements, a DC bias volt- 
age, Vb, generates current, Imw, via resonant tunneling 
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through the QD. (Only a single energy level is relevant 
because the energetic difference between neighboring QD 
energy levels in this device is sufficiently greater than 
kftT and eVg^.) We assume that electrons traversing the 
QD do so ballistically and elastically. As such, electrons 
exit the QD into the down-current electron reservoir at 
an energy higher than the down-current electrochemical 
potential^ Exiting electrons then thermalize within an 
inelastic mean free path, resulting in Joule heat^ This 
effect is added to our FEM model by uniformly depositing 
heat, = ^nw-^Q d ' the electron system within 

one electron-electron mean free path 55 of the QD. Here, 
Rqd is the calculated electrical resistance of the QD, 
which depends on the measured resistance of the HNW, 
and C e . 

We model the phonons in the QD by calculating an 
effective phononic thermal conductivity based on the De- 
bye model^ and the geometry of the InP/InAs/InP het- 
erostructure. The acoustic-mismatch model 3 -^ predicts a 
phonon transmission coefficient of 0.999 at each of the 
InP/InAs interfaces, such that phonon thermal bound- 
ary resistances are negligible at these interfaces. 

The optimization procedure we lay out in Sec. [V] re- 
quires numerous simulation runs with numerous sets of 
HNW parameters. Therefore, instead of modeling the 
Si-SiC>2 substrate, Au leads, and HNW all at once, we 
reduce the computation time by splitting the FEM sim- 
ulations into two separate sub-models: one for the Si- 
SiC>2 substrate and Au leads, and one for the HNW. 
The substrate simulations were ran once and used as the 
boundary conditions for T c and T p h in the HNW sim- 
ulations. As a check, a representative set of HNW pa- 
rameters were simulated both in the HNW model and a 
comprehensive model that includes the substrate, leads, 
and HNW. A comparison of these two models shows a 
difference of roughly 1%, which validates the use of the 
two sub-models for the system. Schematics of the two 
systems with their respective boundary conditions are 
shown in Fig. [3j 

We estimate that radiative losses, and direct heat ex- 
change between the leads/HNW and helium (He-3) bath 
can be neglected when compared to diffusive heat flows 
and the heat exchanged between the electron and phonon 
systems. We also note that ballistic electron and phonon 
effects in the HNW would invalidate Eqs. (TTJ) and @; 
however, these effects are beyond the scope of our paper. 



IV. HEATING THE DRAIN ELECTRONS 

Our model contains three heat flow mechanisms 
through the HNW system capable of delivering heat to 
the drain electron reservoir: (i) Joule heat and electron 
diffusion through the HNW, (ii) phonon diffusion be- 
tween the HC and CC through the substrate, and (hi) 
phonon-mediated heat flow through the HNW. As we will 
show in the following, only mechanism (iii) can account 
for the observed electron temperatures. 



Our simulations predict that both mechanisms (i) and 
(ii) result in AT 0j( j ~ 1 mK, 2-3 orders of magnitude 
smaller than seen in our experiments. The magnitude of 
mechanism (i) is small due to the high electronic ther- 
mal resistance of the quantum dot, which prevents heat 
from reaching the drain electron reservoir. We should 
note that we expect the Wiedemann-Franz law (WFL) 
to break down in our Q D 24 ' 38 ' 39 ; however, the effect of 
mechanism (i) is so small that we can safely ignore Lorenz 
numbers less than lOOLo- The small magnitude of heat- 
ing due to mechanism (ii) is a result of the low heating 
power used. Systems where mechanism (ii) is utilized to 
generate thermal gradients typically have thicker silicon 
dioxide layers and generate significantly more (at least 
two orders of magnitude) Joule heat per unit volume in 
the heating contact^ 

Our simulations predict that mechanism (iii) is capa- 
ble of heating the drain electron reservoir by 100s of mK. 
This mechanism works in three steps to bypass the elec- 
trically insulating QD: (1) transfer heat from the hot 
source electron reservoir to the source phonons though 
e-ph interaction, (2) heat the drain reservoir phonons 
through phononic heat diffusion, and finally (3) use the 
warm drain phonons to heat the drain electron reservoir 
through c-ph interaction. A similar type of energy trans- 
fer process has been previously seen in 2DEG devices^ 

The effects of each heat flow are illustrated by mod- 
eled electron and phonon temperature profiles along the 
HNW axis; Fig. 2] shows three simulated temperature 
profiles for a drain electron temperature measurement. 
Here, Joule heating occurs to the left of the QD because 
for a drain-electron temperature measurement, electrons 
flow from the drain reservoir to the source reservoir i ' 25 
The three temperature profiles shown are the decoupled 
and coupled cases, and a transition state. The decou- 
pled state refers to the state where the electron and 
phonon temperatures are completely independent of each 
other. Likewise, the coupled state refers to when the 
electron and phonon systems act as one effective system, 
T c = Tph- Mechanism (i) results in slight increases in 
AT e ( s d ), and is easily seen in the decoupled tempera- 
ture profile around the QD. Mechanism (ii), seen as an 
increased ATcc, is too small to be visible in the graph 
and has a negligible effect on AT Cj d. Mechanism (iii) 
leads to a transfer of heat between the source and drain 
electron reservoirs, the results of which are the transition 
and coupled temperature profiles. 



A. Two Degrees of Freedom 

To gain a deeper understanding of the factors that de- 
termine the temperature profiles, we now look at results 
from the linearized form of the diffusive heat model de- 
fined by Eqs. (TTJ and @. Importantly, in this model 
there are only two physical degrees of freedom that de- 
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FIG. 4. Example coupled (dot-dash lines) and decoupled 
(solid lines) electron (thick) and phonon (thin) temperature 
profiles through the HNW for a measurement of AT Cj d- Also 
shown is a temperature profile (dashed line) between the cou- 
pled and decoupled regimes, the transition regime. For the 
plots shown here, A ~ 0.9; note that this value was chosen for 
demonstration purposes. The hashed region along the HNW 
position represents the down-current Joule heating due to the 
voltage drop across the QD. The inset shows how AT ei ( Sid ) de- 
pend on the degree of coupling. 



scribe how heat flows through the material: 



^c-ph 



^c-ph 
K ph 



where 



and 



^e-ph 



TnT^d 2 



A = 



K ph 



(8) 



(9) 



(10) 



Here d is a characteristic length scale, see next section, 
and K e -ph can be interpreted as an effective thermal con- 
ductivity between the electron and phonon systems due 
to e-ph interaction. Note that both 8 and A arc di- 
mensionless quantities, and that 6 is derived by solving 
the linearized forms of Eqs. ([1} and ©. 5 represents 
the degree of thermal coupling between the electron and 
phonon systems: 5 <C 1 represents the decoupled state, 
while 8 ^ 1 the coupled state. See the Appendix for 
further details on how 8 and A are used to match the 
simulated AT c ( s( j) to the observed values. 

Since we arc interested in matching our simulations to 
two experimental temperatures, T ( s d), $ an d A are the 
two parameters that we ultimately want to manipulate. 
Due to this limited number of relevant parameters, mul- 
tiple combinations of the original four parameters, {C c , 
C p h, I\ n}, exist that will result in the same 8 and A pair. 
As such, we cannot expect to extract unique values for 
all four physical parameters. 



V. OPTIMIZED PARAMETER SETS 

We search for optimized parameter sets that mini- 
mize the error between the simulated electron temper- 
atures and the experimental data by repeatedly apply- 
ing the Nelder-Mead optimization method^. Each time 
the method is run, the starting Nelder-Mead simplex is 
initialized with random values chosen from within the 
parameter bounds defined in Table [Til but is then free 
to explore an unbounded parameter space for the dura- 
tion of the optimization run. We define the simulation 
error as the RMS difference between the simulated and 
measured temperatures, 




T 



(i) 



c, cxp. 



(11) 



where i denotes the z'th experimental heating current, 
and N is the number of heating currents used in the op- 
timization procedure for each To. We use N = 4 and fit 
to the In < 250 /iA portions of the experimental curves in 
Fig- HI except for To = 4.25 K where we use all four data 
points. After finding numerous parameter sets that fit 
our data for each To, we find two relevant trends between 
the four parameters: one between the conductivity coeffi- 
cients, {C e , Cph}, and one between the two e-ph coupling 
parameters, {n, T}. We also find interesting trends in the 
calculated 8 and A values. These trends are discussed in 
the following. 



A. Optimized Conductivity Parameters 

Here we consider two calculated values from the op- 
timized parameters sets: the ratio of the two thermal 
conductivity parameters, C c and C p h, and the calculated 
A values, the former is plotted in Fig. [5] In both cases 
we find a temperature dependence that scales as a power 
law in Tq~. 



^- = (9.55 x 10 6 

Cph 



W _1 fi _1 K 



-o.i 



)T 4 - 



and 



A = (0.113 KT 2 ' 7 )T 2 - 7 



(12) 



(13) 



Equation (|T2j) demonstrates the presence of an 
unaccounted-for temperature dependence in our model, 
since C' c and C p h were originally assumed to be temper- 
ature independent. Without further measurements how- 
ever, it is unclear whether this temperature dependence 
comes from the electron system, the phonon system, or 
a combination of the two. 

Equation (|13[) implies that the electron and phonon 
contributions to the total thermal conductivity in the 
HNW are comparable between 1.2 and 4.2 K. Addition- 
ally, the electronic contribution is the larger of the two 
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FIG. 5. Scatter plot of the optimized C e and C p h parameters 
with linear fits to each To parameter set. 



for most of this temperature range, with a crossover tem- 
perature of roughly 2.2 K. 

One possibility that may help explain these results is 
a breakdown of the WFL in the InAs portion of the 
HNW. We can check the WFL for our HNW by insert- 
ing Eqs. (HU and ([13} into Eq. (JTUJ) and solving for the 
Lorenz number, L: 



A 



a 



L = L (0.486 K _u - b )T u 



0.6 



(14) 



The apparent temperature dependence of L suggests that 
the WFL does indeed break down in the InAs portion 
of the HNW. Since heat transfer between electrons and 
phonons is the key component to transferring heat be- 
tween the two electron reservoirs in our HNW, it may 
follow that inelastic e-ph scattering is the cause of the 
breakdown; however, the breakdown could also be due 
to other effects, such as changes in the Fermi energy. 

The last point we note about the thermal conductiv- 
ity parameters is that the slopes of the four data sets in 
Fig. [5] are dependent on /(T e ), L, and n p h. If we change 
any of these parameters to remove the observed temper- 
ature dependence in the slopes, only the coefficients in 
Eqs. ([P2"]) -(fl"4" |) change; the general results of our paper 
remain the same. 



B. Optimized Electron-Phonon Parameters 

The relation between n and T resulting from our model 
is shown in Fig.^a). Our results show a band of values, 
from which no particular pair of n and T can be specified. 
Additionally, we cannot use the simulation error to pick 
out any regions along the band since no minimum exists. 

Though we cannot extract unique values for n and T, 
we can estimate the degree of coupling, 5 by inserting 
each optimized parameter set into Eq. ((SJ. Since we are 
interested in e-ph coupling over the length scale of the 
InAs portions of the HNW, we chose d to be half the 
length of the HNW, d = 630 nm. Figure [BJb) shows the 
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FIG. 6. (a) Comparison of the optimized n and F parameter 
sets from the FEM simulations, theoretical calculations, and 
literature valuea^-^. The F values here have been normal- 
ized by the carrier density to allow for a comparison similar to 
the power loss per carrier. The theoretical data shown here 
include both deformation and piezoelectric e-ph scattering; 
each theory data point represents a different electron temper- 
ature between 0.7 and 4.7 K and a background temperature 
of 2.94 K. (b) Scatter plot of the <5 values, calculated from the 
optimized FEM parameter sets, as a function of To. 



calculated 6 values as a function of To. The range of 
5 values implies that the electrons and phonons in the 
HNW are essentially in the fully coupled regime. 



VI. THEORY 

We compare the optimized parameter sets extracted 
from our simulations with theoretical equilibrium power 
loss calculations in a nanowire. As a model, we use a 
cylindrical nanowire of radius R = 25 nm with an infi- 
nite potential well in the radial direction, assuming plane 
waves with a parabolic dispersion along the nanowire 
axis: 



R vj {r) 



1 



Akz 



^RJ u+ i(a,vj) ' 



(15) 



(16) 



where a V j is the j-th zero of the Bcssel functions J v {x). 
In order to determine the chemical potential and thus 
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the occupied electronic sub-bands, the three-dimensional 
doping density is translated into the one-dimensional 
electronic density: nid = n^d^R 2 . With this, we can 
calculate the average power loss rate per electron^ 3 - 



P{T e 



At N ^ 



dN q 
dt 



(17) 



where hujq is the phonon energy, iV q is the occupation of 
the phonon mode q, and N is the number of electrons 
in the normalization volume considered. The transition 
rates dN^/dt are evaluated by Fermi's golden rule for the 
scattering mechanisms addressed below. 

For the interaction of the electronic system with the 
lattice, that is, the e-ph interaction, 



v'j'k' ujk 



v 3 k 9v'j',vj 

q 



(fet+6J, (18) 



we consider confined one-dimensional phonons corre- 
sponding to the modes present in an isolated nanowirc. 
These are calculated within an isotropic continuum 
model following Ref.|43|and0, and assuming free-surface 
boundary conditions. This results in a multitude of quan- 
tized phonon modes w qK with k labeling the individual 
modes. 

We consider only comprcssional modes in the calcula- 
tions; for the deformation potential coupling, torsional 
modes do not couple. Flexural modes in principle cou- 
ple to higher, azimuthally dependent states, and are ex- 
pected to increase the power loss beyond the values cal- 
culated here. For piezoelectric coupling, the piezoelec- 
tric constants for wurtzite InAs are not known. We thus 
calculate the constants from the zinc-blende value, eu, 
using the transformation described in Ref. l45l . There, a 
zinc-blende lattice in the (111) growth direction is consid- 
ered, which corresponds to the wurtzite lattice in a near- 
est neighbor approximation. The parameters used for our 
calculations are the lattice temperature To = 2.94 K, the 
deformation potential of the conduction band D = —5.08 
cV, the mass density p = 5680 kg/m 3 , and the longitu- 
dinal speed of sound vi = 4410 m/s. For more details, 
see Ref. H^. We note that this theory predicts equal con- 
tributions from the deformation and piezoelectric e-ph 
couplings. In bulk materials, optical phonons are typ- 
ically assumed to contribute negligible amounts due to 
being frozen out. 

Using Eq. (|17[) in combination with Eq. we can 
convert P(T C ) to an nvsT plot, 



n = T, 



dP/dT c 
P(T C ) 



T"n, T^n 

1 c J ph 



(19) 
(20) 



The result of these calculations are shown in Fig. Eta) 
alongside our FEM results. 



Interestingly, the shape of the curve mapped out by the 
theory values closely follows the shape of our FEM re- 
sults instead of localizing about a single (n, T) pair. The 
range of n and T values in both data sets, particularly 
the theory data set, suggests that e-ph interaction in our 
HNW does not follow the power law equation given by 
Eq. ([3]); previous research on similar NWs found an expo- 
nential dependence 1 ^. Also of note is that the theory val- 
ues are around 10 3 times smaller than the FEM results. 
This observation, that deformation potential scattering 
with acoustic phonons underestimates the energy trans- 
fer, goes well with previous transport studies in gated 
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nanowire&aii. 



VII. COMPARISON WITH LITERATURE 
VALUES 

Finally, we compare our results to the electron energy 
power loss experiments conducted by Sugaya et alpii and 
Prasad et a h 12 ' 13 Note that the wires in both of these 
papers were constructed out of InGaAs 2DEG materials, 
compared to our vertically grown HNW. The 2DEG wires 
in these papers range from 25 to 770 nm wide and 8 to 
25 nm thicke r 13 ' 47 ' 48 The n and T/n carr i er values from 
their data are shown in Fig. [f^a). 

To within a factor of 10, these values agree well with 
our parameter sets, with our band of values being the 
larger. Theoretical calculations by Vartanian et ali^S 
have recently shown that increasing the degree of lateral 
confinement in NWs enhances the electron energy loss 
rate. Since our HNW has a larger degree of confinement, 
due to being vertically grown instead of etched out of a 
2DEG, it should then exhibit larger electron energy loss 
rates, which is exactly what we see in Fig. [6ja) . 

It is interesting to note that, like our theory values, 
the literature values lie along a band similar in structure 
to our FEM results. Whether this is indicative of an 
underlying physical trend, or merely a consequence of 
the equational form of Eq. ([3]) is beyond the scope of this 
work. 



VIII. DISCUSSION AND CONCLUSION 

We have shown that by combining FEM with exper- 
imental temperature measurements, we can investigate 
the thermal properties of InAs nanowires, including e-ph 
interaction. We find that the polynomial form for e-ph 
interaction given by Eq. ([3]) doesn't match our theory 
calculations, earlier experimental results, nor our data. 
We have also demonstrated that between 1.2 and 4.2 K, 
the electrons and phonons appear to be in a coupled state 
with comparable thermal conductivities. 

The large difference between the FEM and theory data 
sets in Fig.^a) demonstrates that additional thermal ef- 
fects, either in our FEM model or in our theoretical calcu- 
lations, are missing. This is true even with the theoretical 
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result of enhanced e-ph scattering rates due to piezoelec- 
tric scattering, as observed in previous experiments^. 
In our FEM model, we have neglected thermal boundary 
resistances, electron-phonon Kapitza conductances^— , 
and non-equilibrated electron and phonon distributions. 
Due to our system being in a relatively coupled state, 
these additional thermal effects could affect the extracted 
5 and A values, and thus the (n, T) band in Fig. IHJa), re- 
ported here. From the theory side, we previously noted 
that the inclusion of flexural phonon modes in the HNW 
could increase the predicted coupling strength. It is un- 
likely that any single effect will reconcile the roughly 
three orders of magnitude difference between the theory 
and simulation data in Fig. [6ja). Instead, a combina- 
tion of additional thermal effects in the FEM model and 
refining the theory calculations is likely required. 
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APPENDIX 

Here we look in greater detail at how 5 and A affect 
the electron and phonon systems, as well as how they 
can be used to match the simulated AT e ( s d ) values with 
experiments. 

If we interpret the Kj's in Eq. © as heat fluxes instead 
of conductivities, 5 can be thought of as a comparison of 
the rate of heat exchange through e-ph interaction, to 
the rates at which the electron and phonon systems in- 



dividually dissipate heat to their surroundings through 
diffusion. If cither system dissipates heat slower than 
the rate of e-ph heat exchange, K e -ph/K( CiP h) 3> 1, then 
that system can be considered coupled. For example, if 
one has K -ph/^c 3> 1 and K e -ph/^ph "C 1, T e will adjust 
to match T p h- Meanwhile, T p h will remain unchanged 
from its decoupled state since the phonon system can dis- 
tribute heat much faster than the e-ph energy exchange 
rate. As this example demonstrates, one could conceiv- 
ably have the case where one system's temperature be- 
haves as though it is coupled, while the other behaves as 
though it is decoupled. 

Alternatively, if both the electron and phonon sys- 
tem can be considered coupled, K _ p h/fto ^ 1 and 
Ke-ph/ftph ^ 1; A can be used to determine which, if 
cither, of the two systems will behave in an uncoupled 
manner. For example, consider a point along a wire 
where A < 1, n c - P h/^c S> 1, and «e-ph/«ph > 1. If 
T c > T p h, say due to Joule heating, then the net effect of 
c-ph interaction will be for the electron system to rapidly 
transfer heat to the phonon system. Since k c <C K p h, 
the phonon system will be able to dissipate its gained 
heat to the surrounding material at a faster rate than 
the electron system can replenish its lost heat from any 
heat sources. A similar argument applies to the case for 
A> 1, where T e will remain unchanged from its decou- 
pled state, and T p h will adjust to match T e . For A values 
near unity, in a homogeneous material the electron and 
phonon temperatures couple to an effective temperature, 
T e ff, that lies between the two uncoupled temperatures: 

^(e,ph) 1 5=0 < T c ff < T( phiO )|5 =0 . 

For our HNW system, in combination with the 
four externally controlled thermal boundary conditions: 
{AT ei Hc, AT ejC c, AT ph! H C , AT p h,cc}, <5 and A determine 
the balance of heat flowing between the source and drain 
electron reservoirs via the phonon system with the heat 
gained and lost through electron diffusion to the Au 
leads. Properly tuned, the two parameters allow the 
simulated electron temperatures around the QD, AT e>s 
and AT Ci d, to be matched to the experimental measure- 
ments. More specifically, A effectively determines the av- 
erage temperature of the QD in the fully coupled state, 
Tqu = (T CtS + T ei <j)/2, while S determines the tempera- 
ture drop across the QD, AT QD = AT 0iS - AT c , d . Both 
of these effects can be seen in Fig. @] and its inset. 
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